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1. Introduction 

Since the possibility that we are living on a brane in a higher dimensional spacetime 
was proposed, much effort has been devoted to investigating the cosmology of brane 
worlds pp. A model proposed by Randall and Sundrum (RS) provided a very simple 
and interesting playground in which to investigate brane world cosmology [2]. In their 
model, our world is realized on a four- dimensional brane in five-dimensional spacetime. 
The remarkable feature of this model is that the size of extra dimension could be infinite. 
A homogeneous and isotropic cosmological solution was found jSJ-jH] and then the theory 
of cosmological perturbation has been extensively investigated [Zj-|SZ| (see [SB] for a 
review). The study of cosmological perturbation is very important because it can 
provide a means for testing brane world models using forthcoming precise cosmological 
observations. However, it turns out that the analysis of cosmological perturbation 
in brane worlds is extremely difficult. Although a large number of papers has been 
published on this subject, there are rather few quantitative predications. This is because 
the evolution of perturbations on the brane is inevitably coupled to the perturbations 
in the five-dimensional bulk. Thus we need to solve very complicated coupled partial 
differential equations with complicated boundary conditions that arise from the brane. 

Although some progress has been made in a model where a higher dimensional 
spacetime is bounded by two branes [39J 40J, there is still no quantitative prediction 
of the evolution of perturbations in a single-brane model. The crucial difficulty of the 
single brane model is that the bulk spacetime extends to an infinity. A useful method 
for tackling this problem was proposed in Ref. [T2j, (see also [TH], [H] and [HE])- One 
point is that the background 5-dimensional spacetime is just the anti de Sitter (AdS) 
spacetime (or AdS-Schwarzshild spacetime) in the RS model. Thus we can easily find 
general solutions for the perturbations throughout the bulk. A difficulty arises from the 
existence of the brane. We need to select a particular solution that correctly satisfies the 
boundary conditions at the brane. In this paper, we present a solution for this problem 
using a numerical method. Then the late time behavior of the perturbations is derived 
by solving the 5-dimensional perturbations throughout the infinite bulk. In this paper 
the bulk spacetime is assumed to be AdS spacetime without black holes and we only 
consider the late time behavior of the perturbations. The generalization of the analysis 
will be discussed in the conclusion section. 

The structure of the paper is as follows. In Section 2, we describe the set-up of 
our model. In section 3, the solution for background spacetime is derived. In section 4, 
an evolution of tensor perturbation is considered. We describe the numerical method 
used to solve the moving boundary conditions here. In section 5, a more complicated 
evolution of scalar perturbations is derived. In section 6, we try to interpret our result 
using the gradient expansion method. Conclusions and a discussion of the generalization 
of our work are given in Section 7. 
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2. The setup 

We consider the 5D action of the RS model: 

S =J^J d 5 xV=g~ (n 5 + 12/i 2 ) - A J d 4 xV^g~^ e 

-\- J d X\J (jbr-ane^ -matter j (1) 

where TZ 5 is the 5D Ricci scalar, p is the curvature scale of the AdS spacetime and 
k 2 = 87rG 5 where G5 is the Newton's constant in the 5D spacetime. The brane has 
tension A and the induced metric on the brane is denoted as gbrane- The tension A of 
the brane is taken as k 2 X = 6/i to ensure that the brane becomes Minkowski spacetime 
if there is no matter on the brane. Matter is confined to the 4D brane world and is 
described by the Lagrangian C ma tter- We will assume Z 2 symmetry across the brane. 
On the brane, the metric is given by 

ds\ Tane = —dt 2 + a (t) 2 Sijdx l dx j . (2) 

We will decompose the perturbations into scalar, vector and tensor ones defined with 
respect to 3-space 5ij and only consider tensors and scalars in this paper. The energy 
momentum tensor of the matter confined to the brane is given by Tjf = T^ u + 5Tjf\ 

Tbu = diag(—p,p,p,p), p = wp, (3) 

m = f :// ) • (4) 

Here the matter anisotropic stress is neglected for simplicity. 

3. Background spacetime 

In this section, we briefly consider the background spacetime. In this paper, we only 
consider the maximally symmetric bulk spacetime, that is, we assume that the bulk 
spacetime is given by Anti de Sitter (AdS) spacetime without a black hole mass. 
The simplest coordinate system for describing AdS spacetime is given by the Poicare 
coordinate: 

ds 2 = (J-^J (dz 2 + r] liv dx^dx v ). (5) 

Our brane is moving in this coordinate system. The motion of the brane is determined 
by imposing the junction condition. The brane motion is described by [1] 

' =T(t), (6) 



fia (t) 



where 



\a J 3 36 
T =-Jl+(-V, (7) 
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where t is the cosmic time on the brane and a dot denotes the derivative with respect to 
t. The expansion of the universe can be understood as the motion of the brane through 
the bulk. 

Of course, the statement that our brane is moving is a coordinate dependent 
statement. It is possible to choose coordinates where the position of the brane is fixed. 
The Gaussian-Normal (GN) coordinate is such a coordinate. The metric is given by [3] 



ds 2 = dy 2 - n(y, tfdt 2 + a(y, t) 2 Sijdx l dx j , (8) 



where 



a(y,t) = a a (t) 



cosh ay — 1 H sinn ay 

V 6/i I 



( k2 P \ 

n(y, t) = cosh ay — 1 (2 + 3w) sinh ay. (9) 

\ 6fi J 

The brane is located at y — 0. This coordinate is convenient for imposing the junction 
conditions because the position of the brane is fixed. The junction conditions in the 
background spacetime are given by 





(10) 

where the subscript b is used to denote bulk quantities evaluated at the brane. Using 
the Freidmann equation, the former equation can be written as 




\ 



■ (?) 

Note that there is a coordinate singularity in this coordinate y c at the finite distance 
from the brane where a(y c ,t) = 0. This hyper surface corresponds to the past Cauchy 
horizon of the AdS spacetime. Thus, this coordinate is not suitable for discussing the 
global structure of the solution in the bulk. However, in order to impose the junction 
conditions, we need only know the geometry of the bulk near the brane. Thus we will 
solve for the bulk perturbations in Poincare coordinate and move to the GN coordinate 
when we impose the boundary conditions. In order to do so, we need to rewrite the 
junction conditions in the GN coordinate as conditions in Poincare coordinate. The 
coordinate transformation between the two coordinate systems has been investigated in 
Ref.[H]. The transformation is very complicated, but we only need the formula relating 
the derivative with respect to the GN coordinate to the derivative with respect to 
Poincare coordinate on the brane. This is given by the following formulae; 
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Let us begin with the tensor perturbation: 



ds z 



dr + (5ij + h i j)dx l dx :i ) . 



We expand hi* as 



hij 



d 3 k 

J2^y 



h(z, T)eije ikx , 



where is the polarisation tensor. The evolution equation is simply given by 



d 2 h 3 dh d 2 h 



dz 2 z dz dr 2 
It is easy to find a solution 



- k 2 h = 0. 



(13) 
(14) 

(15) 

(16) 
where 

Z 2 (mz) = Hf\mz) + b{m)H^\mz), (17) 

and to 2 = m 2 + k 2 . Here and are Hankel functions of the first kind and the 
second kind respectively and h(m),b(m) and c(m) are arbitrary coefficients. This is a 
general solution for tensor perturbations in the AdS bulk. Our spacetime is bounded 
by the brane. Thus we should impose a boundary condition via a junction condition. 
In the GN coordinate, the junction condition is simply given by 

'dh' 



h(z,r) = / dmh{m)z 2 Z^iviiz) {e lu)T + c(m)e 



dy 



0. 



;is) 



Using the formulation equations (fT2|) . we can rewrite this condition in terms of the 
Poincare coordinate 




-1/2 



dh 



0. 



Z=l/fMo(t),T=T(t) 

Let us impose the boundary condition Eq. (jl9J) on a general solution given by 

dmh(m)m 2 z 2 Z2(mz)e luJT . 



h(z,r) 

The boundary condition of equation (flUjl gives 

,3 



dmh{m) 



0. 



m 
H 2 a 2 



m 



jja Q 




-1/2 



lujm 
y?a 2 Q 



m 



fia 



(19) 



(20) 



(21) 



The problem is finding the solution for h(m) that satisfies this boundary condition. In 
order to find the solution, we need a numerical method. 

As suggested by recent numerical works [12], [3*5] . we expect the behavior of the 
perturbation to be well described by a standard 4D Einstein gravity at low energies 
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H/ /i <C 1. We can argue on the recovery of the 4D solution as follows. At late times, it 
is natural to assume that the Kaluza Klein (KK) mass is small compared with the bulk 
curvature scale because the Hubble scale of the brane universe is much smaller than the 
bulk curvature scale H/ /i < 1. So we assume 

— = e « 1. (22) 
fia 

Then expanding in terms of e by using the asymptotic behavior of Hankel functions 
with small arguments, 

H , H ^Jj^Jl(lY, mz) ^t^l(l)\ (23) 

7T \Z J 7T \Z J 

the general solution of equation (j2*U|) on the brane becomes 

h b (t) ~ J dmh{m)e iujri , (24) 

where we neglected the numerical factor and used T(t) = 77 at late times where rj is the 
conformal time. On the other hand, the boundary condition of equation (|21|) becomes 



Q = J dm -2iuH + (e 4 ,e 4 log e)\ 



h(m)e 1 ^. (25) 



Then on rewriting m 2 = k 2 + uj 2 and uj in terms of the derivative with respect to 
conformal time, this condition can be written as 

( d 2 1 da Q d , 9 \ , . . 

-2 — ^- k ) h b = 0. (26) 



\ drf a Q dt] dr] 

This is noting but the evolution equation obtained in 4D Einstein gravity for the tensor 
perturbation. It implies that, at late times H/fi <C 1, the boundary condition selects 
a particular solution that obeys the 4D evolution equation on the brane. It should be 
noted that the excitation of KK modes in Poincare modes is necessary in order to satisfy 
the boundary condition because zero-mode solution is just e lkv . The movement of the 
brane excites KK modes. They give the damping of zero-mode and give the friction 
term in Eq. (j26|) . 

However, the above argument assumes that the junction condition (|21j) can be 
satisfied by the KK modes with the condition (}2*2*j) . It only ensures that once there 
is a solution for h(m) that satisfies the condition (|2*2*|) . the 4D evolution equation is 
recovered. We should check that there exists a solution for the junction condition (|2*T|) 
that satisfies the condition (J22j). For this purpose, the spectrum of the KK mass should 
be determined by imposing the junction condition at the brane (j21)l as well as the initial 
conditions and the boundary condition in the bulk. 

Let us try to find solutions for b(m),c(m) and h(m) numerically. We first specify 
the boundary condition at the AdS past Cauchy horizon, which corresponds to the 
infinity z — * 00. The most natural assumption is that the wave is outgoing, so that 
there is no incoming radiation. We demand 

hoc [ dmz 3/2 e imiz - T \ (27) 
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for z — > oo and k — > 0. This condition gives 6(m) = c(m) = 0. Then we impose the 
boundary condition of equation (j21j) at the brane. This condition is formally written as 



We numerically find solutions for h(m) that satisfy this boundary condition 
approximately. We first discretize time t and a KK mass m. Then we end up a matrix 
equation. A problem is finding eigenvectors with a null eigenvalue. In principle there 
could be an infinite number of solutions for h(m). However we expect that an excited 
KK mass is small at low energies Hj /i C 1 as long as we consider an initial condition 
that does not contain a significant excitation of KK modes. Thus we introduce a cut-off 
for m. Then we can find finite numbers of approximate solutions for h(m). In addition, 
we should determine the initial condition. Here we adopt an ad hoc initial condition 
that the perturbation is constant with respect to t on the brane at super-horizon scales. 
This can be achieved by mixing the real part and the imaginary part of h appropriately. 
Then we obtain a real solution for h. In appendix C, the accuracy of the numerical 
calculation is shown. 

Figure 1 shows two solutions for h(m) from numerical results. We take H/fi = 
0.0001 at the horizon crossing as an example. As expected we need the KK modes to 
satisfy the boundary condition but the excited KK mass is small compared with /i. If 
we sum the mode functions with the weight h(m), we recover the 4D behavior of the 
perturbation on the brane (Figure 2). 

We would like to make a comment on the cut-off of the KK modes introduced in the 
numerical calculations. As is seen from the Figure 1, the solution obtained for h(m) is 
localized around m/yU < 0.2 which is well below the artificial cut-off at m/fi = 1.2. We 
have confirmed the existence of the solution for h(m) that is localized around mj fi < 0.2 
even if we increase the cut-off up to m//i = 2.4. Thus the existence of the solution with 
mj ' \ia Q << 1 does not depend on the artificial cut-off of the KK mass. On the other 
hand, the introduction of the cut-off certainly kills the solutions that have large mj [ia Q . 
Thus our analysis is limited to the solutions with the KK mass below the cut-off. 

The different solutions for h(m) correspond to different initial profiles in the bulk 
(Figure 3). In this sense, we have not solved the problem yet. In order to find the 
solution that corresponds to a given initial profile, we again sum solutions with different 
h(m). Fortunately, we find that the different solutions h(m) give the same solution on 
the brane. This indicates that the recovery of the 4D behavior does not depend on 
the choice of the initial profile, as long as the initial profile does not contain large KK 
masses. 

If we include the high energy corrections 0(H/ /i), the behavior of the perturbation 
significantly deviates from 4D one. The difficulty of the high energy corrections is that 
the corrections depend on the choice of the initial conditions. Unlike the late time 
evolution, the different solutions for h(m) give different evolutions on the brane. Thus, 
unless we perform the summation of the solutions with different h(m) we cannot get the 
final answer. As we mentioned, the different solutions for h(m) correspond to different 




(28) 
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Figure 1. Solutions for h(m). The vertical line is m//i. We take H//J, — 0.0001 at the 
horizon crossing and the universe is assumed to be radiation dominated. The cut-off 
is introduced at m/ /i = 1.2. 




Figure 2. Behavior of the solution hb on the brane. The horizontal axis is the 
cosmic time. The solid line is the prediction of 4D cosmology with the same initial 
condition. The points represent the solutions in the brane world obtained by summing 
the Poincare mode functions with the weight h(m). All solution for h{m) found in 
numerical calculations reproduces the 4D behavior of the perturbation. 



initial profiles. Thus this indicates the sensitivity of the solution to the choice of initial 
conditions. A detailed study of the high energy correction is beyond the scope of this 
paper and it will be reported in a separate publication |%T] . 
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Figure 3. The initial profile of the solution h(U,z) in the bulk. The horizontal axis 
is jiz. The solid line corresponds to the solution derived from h(m) in the left panel 
of Figure 1 and the dotted line corresponds to that in the right panel. 



5. Scalar perturbation 



Now let us consider the scalar perturbation. There are several ways to calculate the 
scalar perturbation in AdS spacetime. Here we describe an approach based on a master 
variable. Using the generalized 5D longitudinal gauge, the perturbed spacetime is given 
by HU 

ds 2 = b(y, t) 2 (l + 2A yy )dy 2 + n(y, t)A y dydt - n(y, t) 2 (l + 2A 2 )dt 2 



+ a(y, t) 2 (l + 21Z)6 ij dx i dx j . 



(29) 



It was shown that the solution for metric perturbations can be derived from a master 
variable Q [22] [ZZ]; 



.4 
A 
A 



6a I b 2 



mi 













na 




n 






tt" - ^ 


t\ 




ba 








n" + ( 


t\ 




ba 







\ b n 



+ 



n- 



n-\2- 



n 
n 



n 



-fi 2 n 



n' V 
n b 

n b I 



+ 



n- 
1 



2Q 



-n 



n 



?? 



n 



+ fi 2 n 



n b j 



2/i 2 fi \ ,(30) 



as long as the master variable Q in the bulk satisfies a wave equation given by 



n 



na° 



{'ba'' 



of 



nb 



Q = 



(31) 
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In Poincare coordinate, the master equation is given by 
d 2 Q 3dtt fj? d 2 Q 



+ 737 + ^-^ + V^ = (32) 



dz 2 z dz z 2 dr 2 
The solution is easily found as [H] 



2 r d 3 k 

n = — p~ 3 / j— / dmS(m)Z (mz)e- iuJT e ikx . (33) 



z^ J (2tt) 3 

A factor —2 was added just for later convenience. 

The junction conditions are much more complicated. In order to impose the 
junction conditions, we move to the GN coordinate. The junction conditions have 
been found already in literatures (see [S] and |2Tj). but, for completeness, we briefly 
discuss the way to impose junction conditions on a master variable in Appendix A. The 
junction conditions give the expressions for matter perturbations on the brane in terms 
of a master variable; 

a o ^ Sp = _ *! (tf _ O _ 3 ^ (ft _ VL&\ , 
a 2 \ a J a \ n J 

a Q K 2 5q = — ( Q' f2 ] , 



n 

a' - a I • , n' ■ \ I a l a' n' \ ( a \ ' / n 



2 j 



-I 


/a' 


n' 




\ a 


n 



a K 2 5p = fi' — n + 2- m 7 + < 4- + 2 - - - ^fi 

a a\ n J [a \a n J \aj \nj J 

a 2 \a n J \a n J \ a n J 

(34) 

where the prime denotes d y and dot denotes dt in the GN coordinate and the right hand 
side of the equation should be evaluated at the brane. Imposing the equation of state 
5p = c^Sp, we get the boundary condition for fl. 

It is more convenient to rewrite these equations into the follwoing form using the 
expressions for metric perturbations in terms of a master variable equations (j3*Uj) as 

- l -K 2 H5p = 3H(TZ - HA) + \k\8 P£ , 

-l K 2 H5p = - it - H (3 - ^Ak + HA + (2H + 3H 2 - ^-j A 
lk 2 A 1/ H\k 2 ^ 1/ H\ e 

-3^-3 C 1 -^; ^ n -- Q [ i+ H 2 ) K ^ 

-^K 2 H6q =K-HA--i$6q e , 

=\(A + K) + 45ir £ , (35) 

where we defined 

H=[-\ -,< 



a! 



\ 



1 + (j) , k\ = K 2 p, (36) 
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and 



k A n 

k 2 fd 
3al \a Q 
1 

Qal 



i^ld-Ke = — ( 3Q - 3HQ + — Q - 3 



72' 



n' a' 



n a 



n'\. 



(37) 



These equations can be derived from the projected Einstein equation on the brane 



Guv + £, 



^7Tp,v + ^ n 



where 



n 



--T, a T« + + ^(dTapT* - (T:) 2 ) 9i 



(3f 



(39) 



and is the projected 5D Weyl tensor. Here 5ps,5p£ and 5qe are perturbations of 
"Weyl fluid"; 



5E» = k 2 ( Z 5f>e 6q£ ' 1 

4 V a o 2 S a £,i fiPS $ij + fas . 



(40) 

Substituting the solution for Q (|33|) . The metric perturbations and Weyl tensor are 
obtained as 

2 



K 



dmS(m) 



m i m 



pa 



pa 



1 

+ - 



3 \pa 



k \ i in 



pa 



-iwT 



(41) 



2k 4 
3a 



i ! dmS(m)p 2 Z 



in 



-uoT 



pa a 



Ik 2 r 

k 2 J><1£ = 7^ J dmS(m 

2 2v 

a o K 4 0ii £ = 

- H 2 pT 2 I dmS(m) 



m 



K pa D/ 

-2-2 - 



Hmp- 3 Z 1 



711 



rn 



dmS(m)-(2k + 3m)p a Q Zq 

3 ' \pa 



-iu)T 



in 



pa a 



-Z x 



rn 



ua 



- {k 2 + 2m 2 )p- 2 a- 2 Z 



in 



pa Q 



+ 



2HHp 2 J dmS(m)iujmp 2 a Q 2 Zi 



rn 



pa 0/ 



-iuiT 



It should be noted that these equations have been derived already in Ref . [23] . In Ref . [ 
we solved the perturbations in Poincare coordinate using the Randall-Sundrum gauge 
(see Appendix B). Of course the final result completely agrees. 

Now it is possible to write the equation 5p — c 2 5p = in terms of the soltuion for 
Q. At low energies H/p <C 1, the equaiton is simplified very much because —k 2 T~C = k 2 
and we can neglect the terms proportional to H/Ti 2 . For cosmological problems, it is 
natural to assume that the physical 3D wavelength of perturbations is much larger than 
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the AdS curvature length at low energies. Thus we can neglect the terms suppressed by 
k//j,a . Then we end up the 4D Einstein equation except for the equation 

A + H = -a 2 K\5'Ke, (43) 



where 



1Z = — J dmS{m) 



f m 






)■ 


\fia 






m 


m) ^ 






jj,a 




a z K%5>K£ = / dmS(m) Z e~ lujri . (44) 



Here we implicitly assume that there is no dark radiation perturbation (see the 
conclusions section for a discussion of the dark radiation perturbation). At super-horizon 
scales, the conservation of curvature perturbation on the hypersurface of uniform energy 
density; 



C = K--7r\T;- A )> ( 45 ) 




can be shown without using the equation that relates TZ and A [T7|. However, in order 
to predict the Cosmic Microwave Background (CMB) anisotropy we need the relation 
between TZ and A because the large scale Sachs- Wolfe effect is given by 
AT 

— = C + A-K. (46) 

Thus we need to evaluate a Weyl anisotropic stress Stts in order to address the CMB 
anisotropy. 

As for the tensor perturbations, if we assume 

<=(—)<i. <«) 

\fia J 

and using the asymptotic formula for Hankel function 

Hq 2 \z) ~ — i— (7 + log ~ \ + 1, 7 : Euler number, (48) 
7T V 2/ 

we get TZ = —A Then we recover the 4D Einstein gravity. However, as in the case 
for the tensor perturbation we should determine S(m) to justify this assumption. The 
problem is the same as the tensor perturbation one. We used the same boundary and 
initial conditions. The cut-off of the KK mass is introduced in the numerical calculations. 
As for the tensor perturbations, the obtained soltuion for S(m) is localized well below the 
cut-off. In addition, for scalar perturbations, we neglect the terms suppressed by H/jj, 
and k/fia Q and use the 4D Friedmann equation in performing the numerical calculations. 
Figure 4 shows the solution for S(m). From these solutions we can construct the solution 
for metric perturbations. We recover the 4D behavior of the perturbations as expected 
(Figure 5). In particular, the metric perturbations satisfy the relation TZ = —A. The 
behavior of metric perturbations does not depend on the solutions for S(m); thus the 
recovery of 4D cosmology does not depend on the initial conditions. 
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Figure 4. We show two of the solutions for S{m). The horizontal axis is m/fi. The 
perturbations enter the horizon when H//J, = 0.0001. The universe is assumed to be 
radiation dominated. The initial condition is taken such that metric perturbation is 
constant at super horizon scales. The cut-off is introduced at mj [i = 1.2. 




Figure 5. The evolution of the metric perturbation 1Z and —A on the brane. The 
horizontal axis is the cosmic time. The solid line is the prediction of 4D cosmology with 
the same initial condition and the points represent the solutions obtained by summing 
the Poincare mode functions with the weight S(m). The amplitude is normalized so 
that 1Z = 1 at the initial time. 

6. The GN coordinate view 

In this section, we reconsider the evolution of scalar perturbation in the GN coordinate. 
As mentioned in section 3, the GN coordinate has a crucial disadvantage because it 
cannot cover the whole bulk spacetime due to the coordinate singularity. However it is 
still a good coordinate near the brane and we can understand what is happening near 
the brane in a very simple way. Thus it would be instructive to compare the approach 
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used in this paper with the analysis in the GN coordinate. 

Let us consider the low energy/near brane approximation [33] [34] . The metric is 
simply given by 

ds 2 = dy 2 + e~ 2m {-dt 2 + a% 3 dx l dxi). (49) 

Note that this metric is valid only near the brane 

e 2m < (K 2 p/i2)-\ (50) 

At late times, the derivative with respect to time is expected to be much weaker than 
the derivative with respect to y: 

n < n" ~ p- 2 n. (51) 

With this approximation the master equation becomes 

Q" + 2/iJ]' + p 2 Q = 0. (52) 
Then the solution is given by 

n = VL (x)pye-w + C(x) e -' iy . (53) 

Here we give a weak dependence on brane coordinates to the constants of integration 
flo(x) and C{x). These dependences can be determined by imposing the boundary 
conditions. The junction condition is greatly simplified due to the fact a' /a = n'/n at 
late times K 2 p/p <C 1. The junction conditions of equations (|3*l)l are written as 

k 2 

a K 2 5p = -T — ?>HT, 

aj 

a Q K 2 5p = F + 2HT + ET ', (54) 
where we defined 

T = Q' + pfl. (55) 
On imposing the equation of state 5p = c 2 5p, the equation for T is obtained as 

T + (2 + ZcDET + HF + c 2 s k 2 ar 2 T = 0. (56) 
A point here is that T is written only in terms of Qq(x), 

T = /iO (x) (57) 

Thus we cannot determine the function C(x). It is a natural consequence. We are 
basically solving the second order differential equation with respect to y. The boundary 
condition at the brane is insufficient for determining the solution. We should specify 
the boundary condition in the bulk to determine the solution completely. However, the 
GN coordinate is not suitable for this purpose due to the existence of the coordinate 
singularity. Moreover, the gradient expansion method for solving the perturbation 
cannot work for large y (see equation (j50jl ). However, if we calculate the behavior 
of metric perturbations, the leading order term that comes from C disappears; 

A = ^-p 2 n (x) + o(C), 

za, 

TZ= -— p 2 Q {x) + O{C). (58) 
2q 
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Thus if the condition 

0{C) < f/n (59) 

is satisfied, we can ignore the contribution from C(x). Then the metric perturbations 
are solely given by Qq. Hence the equation for T gives the evolution equation for metric 
perturbations 

A = -K, 

71 + (4 + 3c 2 s )H7l + (2H + 3H 2 + 3c 2 s H 2 )7l + c 2 k 2 a- 2 7l = 0. (60) 

These are nothing but the evolution equation in standard 4D cosmology. 

The crucial disadvantage of the GN coordinate is that there is no way to find the 
solution for C as mentioned above. Hence it is impossible to justify the condition (jEHj) 
that is needed to show the recovery of the 4D evolution equations. And also the function 
C could modify the relation between A and 7Z 

A + TZ = 0(C). (61) 

This equation is essential for calculating the CMB anisotropy. Thus we can say nothing 
about CMB anisotropy. 

Let us compare the above arguments with the analysis using mode functions in the 
Poincare coordinate. At late times, T can be calculated as 

T = iiQ = 2/^ 2 / dmS(m)mZ l [ — ] e - *"". (62) 

j Va"W 

From Eq.(jHJ), we can relate fl to 71 as 

Q = -2a fi~ 2 TZ, (63) 

which agrees with Eq. (j58j) . The difficulty that arises in the GN coordinate was that there 
is no way to determine C. In the numerical calculation done in Poincare coordinate we 
did impose the boundary condition and initial conditions in the bulk. Hence the function 
C was determined. The solution Eq. (j53|) indicates that C is the value of Q on the brane 
Qb = C (note that brane is located at y = 0). Thus C is given in terms of the solution 
in Poincare coordinate as 

Vl b = C = -2fi~ 2 a [ dmS(m)Z [ — ) e"^. (64) 
J \{J>a J 

We can also derive the equation that relates A and 1Z from Eqs. (|4~H|) and (jUJ); 

A + 1Z = -a 2 Kl5n £ , 

,2^.2, _ f j mQ / ro \ f ™_\ v ( ™\ e -iur, 



aiKlSne = / dmS(m) Z e" 1 ^, (65) 

which should be compared with Eq. ()61|) . The leading order behavior of Weyl anisotropic 
stress in mj \ia Q < 1 contains the non-local term 

a 2 Kl5n £ oc / dmS(m) ( —) log ( —) e*^. (66) 



Late time behavior of cosmological perturbations in a single brane model 



16 



Thus the correction o? k\8ii£ describes the 5D corrections that are caused by the 
propagation of perturbations into the bulk. This explains why the GN coordinate 
approach can tell nothing about this correction. We should determine the solution 
for the bulk perturbations completely to address the 5D correction. And also it is now 
understood that the physical meanings of the condition (JHH)) is that the 5D effect that 
comes from 5D Weyl tensor o? k?J>'K£ is negligible when determining the solution for 
metric perturbations 7Z and A. 

A point here is that we can quantitatively check the condition using the 
numerical results done in the Poincare coordinate. Figure 6 shows the behavior of Weyl 
anisotropic stress from numerical solutions. The amplitude of the Weyl anisotropic stress 
should be compared with the amplitude of 1Z in Figure 5. The smallness of the amplitude 
of alKl5n £ compared with the amplitude of 71 indicates that the 5D effect is negligible 
when determining the solution for metric perturbations. This also quantitatively justify 
the condition (|59|) in terms of the solutions in the Poincare coordinate. Thus we can 
show the recovery of the 4D physics at low energies H/ ' \x <C 1. 
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Figure 6. Solutions for Weyl anisotropic stress k\o^8it£. The amplitude is normalized 
so that TZ = 1 at the initial time (see Figure 5). The horizontal axis is cosmic time. 
The solid line is the solution obtained from the solution for S(m) in the left panel in 
Figure 4 and the dot line is that from in the right panel 



7. Conclusion 



We solved for the evolution of cosmological perturbations in a single- RS-brane world 
model. Assuming that the bulk is AdS spacetime, the general solutions in the Poincare 
coordinate were used to construct the solution. This allows us to find the solution 
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throughout the infinite bulk spacetime. The junction conditions at the moving brane 
are imposed numerically. Then we obtained the solution for perturbations by solving the 
5D bulk perturbations. We considered a late time evolution of perturbations. At late 
times H/ /i<l, the behavior of the perturbations agrees well with that in conventional 
4D cosmology. 

Our result indicates that it is very difficult to find brane world corrections in 
CMB anisotropy in a single-RS-brane model. The AdS curvature scale is restricted 
to yU -1 < 0.1mm from table-top experiments of Newton's force law. At high energies 
H/n > 1, the perturbations that are relevant to the CMB physics stay at superhorizon 
scales. At superhorizon scales, the curvature perturbation on hypersurface of uniform 
energy density is conserved even in a brane world model (note that there could 
be a correction from dark radiation perturbation; see below). Thus the curvature 
perturbation evolves in the same way as 4D cosmology. However, in order to predict 
the CMB anisotropy we need to know the anisotropic stress induced by 5D Weyl 
tensor in addition to the curvature perturbation. This could give a 5D effect to the 
CMB anisotropy. Unfortunately, at the decoupling time when the CMB spectrum is 
determined, H/ \i <C 1 is satisfied to extremely high accuracy. Our result indicates that 
this 5D effect is too small to be observed. 

However, there are still several possibilities to observe the brane world corrections. 
The first possibility is provided by the high energy corrections that arise when H/ \i 
becomes large. This correction is particularly important for tensor perturbations 
because there is a possibility of directly proving the evolution of tensor perturbations at 
high energies if we can observe the stochastic background of gravitational waves from 
inflation. A difficulty of the calculations at high energies is that the behavior of the 
perturbations depends on the initial conditions. Thus we need an extra step to impose 
the initial conditions. A detailed numerical study of the high energy correction will be 
presented in a separate publication We should note that there has been a very 

interesting attempt to understand the high energy corrections in terms of AdS/CFT 
correspondence jSEj- The logarithmic corrections log(ra///a ) in Eqs. (j25j) and (jfifi)) can 
be understood as corrections due to a coupling to the CFT. Although this approach 
can treat only mild corrections Hj\i < 1, it is very important to understand the 
corrections analytically. It would be interesting to extend the analysis of Ref . to 
scalar perturbations. 

Another possibility arises if we allow the existence of a black hole (BH) in the 
bulk. In the background spacetime, a bulk BH induces the so-called dark radiation that 
modifies the evolution of the universe. If we consider the perturbation, dark radiation 
also has a perturbation and modifies the evolution of perturbations. An interesting point 
is that this effect could be large even at low energies and it could give very interesting 
features in CMB anisotropy jBH] • Unfortunately, this dark radiation perturbation would 
be a non-normalizable mode in AdS spacetime [13], [II]- However, there is a subtlety 
in the argument of the normalizability. It has been shown that the dark radiation 
perturbation corresponds to putting a small BH in the bulk [4~3j . We cannot treat the 
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effect of the BH perturb at ively near the event horizon even if the BH mass is small. So 
it is difficult to discuss the dark radiation perturbation in AdS spacetime background. 
We should carefully investigate this mode in AdS-Schwarzshild spacetime. Recently a 
master variable has been found in AdS-Schwarzshild spacetime [45J. Thus it would be 
very interesting to extend our analysis to AdS-Schwarzshild spacetime and investigate 
the late time behavior of perturbations. 

Finally, we should specify the initial conditions from inflation. For scalar 
perturbations, analysis has been done by assuming that the contribution from bulk 
perturbations can be neglected However we should carefully examine the validity of 
this assumption. Fortunately, during inflation, the equations are simplified significantly. 
Thus it is possible to analyze the behavior of perturbations analytically to some 
extent j33j. We hope that the primordial fluctuations will also have some brane world 
signatures. 

Quantitative analysis of 5D effects on cosmological perturbations in a single-RS- 
brane world has just begun and much things remain to be done. We hope that our 
study provides the first step toward detailed predictions of cosmological observations in 
this model. 
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Appendix A. Junction condition for 



In the brane world, we should perturb the location of the brane as well as the 
perturbation in the bulk |32|- The 5D longitudinal gauge in GN coordinate completely 
fixes the gauge, thus there is no guarantee that the perturbed brane is located at y — 0. 
Thus we should go to the gauge where the perturbed brane is fixed at y — 0. We call 
this coordinate the brane-GN coordinate. Under a scalar gauge transformation 

t — > t = t + St, y^y = y + 5y, x — > x = x + 8x, (A.l) 

metric perturbations transform as 



.4 
A 
E 



yy 



A 



by' , A y = A y + n8t' 5y, B y = —5x' 



1 



A- 5t by 5t 

n n 

— 8x, B = —8t — 
a 2 



n 

^ ^ a ' r a P 
K = K 5y - -5t, 



Sx, 



where 



9AB 



( -n 2 (l + 2A) 
nA„ 



o?B i 

a 2 [(1 + 211)5^ + 2E tij ] 
a 2 B y ^ 



nA v 
1 + 2A 



yy 



(A.2) 



(A.3) 
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The conditions for the GN coordinate A yy = A y = B y = give 



yy • 



by' = A 

St' = - -A y + —Sy, 
n n z 

Sx' = - \s y . (A.4) 

or 

The junction conditions in brane-GN coordinate are given by 

K 2 



A' b = -(2Sp + 3Sp), 
o 

B' b = tf^Sq, 

E' b = 0. (A.5) 
Using the junction condition for E, we get 

E' b = -5x' b = ^6y b = 0. (A.6) 

a o 

Thus the brane location is not perturbed; that is 5y b = 0. The metric perturbations 
in 5D longitudinal gauge can be regarded as the induced metric perturbations on the 
brane. There is a residual gauge freedom in St on the brane. Using this gauge freedom 
we can set 

St b = 0. (A.7) 

Then it is possible to write A' b , R' b and B' b in terms of the metric perturbations in the 
5D longitudinal gauge. Substituting the solutions for metric perturbations in terms of a 
master variable, we can express the matter perturbations in terms of a master variable. 

Appendix B. Scalar perturbation in Randall-Sundrum gauge 

An alternative way to solve scalar perturbations is to use the Randall-Sundrum gauge 
in the Poincare coordinate. This was done in Ref.|T2], [2S1- We can start with the 
perturbed AdS spacetime in Poincare coordinate; 

ds 2 = (J^\ (dz 2 - (1 + 2(f))dr 2 + 2b i dx i dr + ((1 - 2#)<% + 2£ ij ) dx i dx j 

(B.l) 

where 4>,b,^ and E are given by 

h = (fiz) 2 [ [ dm h(m)Z 2 (mz)e- iuJT e ikx , {h = <p, b, E). (B.2) 



., (2vr) 3 . 

Here we used the transverse traceless gauge conditions 
6-3V + V 2 E = 0, 
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(B-3) 



Thus the coefficients h(m) satisfy 

2A; 4 



<p[m) 
b{m) 

E(m) 



3m 2 



..Vk 2 + m 2 k 2 ^ 2 
— 4« — Elm) 



fcV 



3m 2 
j 

-E(m), 



2k 2 + 3m 5 
3m 2 



^ 2 E(m) 



(B.4) 



where E(m) is the arbitrary coefficient. As in the case for 5D longitudinal gauge, it 
is possible to relate these solutions to matter perturbations on the brane by imposing 
the junction conditions. The final result completely agrees with the one derived using 
a master variable with the identification E(m) = S(m). 

Appendix C. Numerical accuracy 

In this section we show the accuracy of the numerical calculations. We have checked the 
accuracy of the junction condition. For the tensor perturbation we evaluate equation 
(|2*Tj) using the solution for h(m) (Figure CI). For scalar perturbation we evaluate 
K 2 /x(5p — c 2 5p) using the solution for S(m) (Figure C2). In the figures, we make the 
equations dimensionless using //. 
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Figure CI. The evaluation of Ea. pijl multiplied /i 1 . The left (right) panel shows 
the result for h(m) in the left (right) panel of Figure 1. 

We need appropriate references to compare with the numerical violation of the 
junction condition. We define the "normalized" error of the junction condition by 
dividing equation (J2T|) by the representative term in the junction condiiton. For the 
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Figure C2. The evaluation of K 2 fi(5p — c 2 5p) multiplied by fi 2 . The left (right) 
panel shows the result for S(m) in the left (right) panel of Figure 4. 



tensor perturbation we divided the junction condition by the first term of equation (|21|) 
(Figure C3). For the scalar perturbation, we divided the equation K 2 fi(5p — c 2 5p) by 



/ dmS{m)u] 2 — -Z\ 
J \ia A 



m 



-iuiT 



(C.l) 



which has the typical amplitude of the terms in the equation K 2 fi(5p — c 2 5p). The result 
is shown in Figure C4. By construction, the normalized error is dimensionless and it 
measures the violation of the junction conditions. Note that the appearance of peaks of 
large errors is caused by the fact that the denominator becomes 0; thus it is an artifact 
of the definition of the normalized error. 
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Figure C3. The normalized error for the tensor perturbation. The left (right) panel 
shows the result for h(m) in the left (right) panel of Figure 1. 



For comparison, the deviation from the 4D solution with the same initial conditions 
is shown in Figure C5 for the tensor perturbation. Within the accuracy of the numerical 
calculations, the solution obtained cannot be distinguished from the 4D solution. 
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Figure C4. The normalized error for the scalar perturbation. The left (right) panel 
shows the result for S(m) in the left (right) panel of Figure 4. 




Figure C5. The deviation from the 4D solution with the same initial condition. The 
solid line corresponds to the solution derived from h(m) in the left panel of Figure 1 
and the dotted line corresponds to that in the right panel. 
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